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Abstract. We investigate the attractive Hubbard model in infinite spatial dimensions at 
quarter filling. By combining dynamical mean-field theory with continuous-time quantum Monte 
Carlo simulations in the Nambu formalism, we directly deal with the superfluid phase in the 
population imbalanced system. We discuss the low energy properties in the polarized superfluid 



j— h state and the pseudogap behavior in the vicinity of the critical temperature. 

a 

'"O ' 1. Introduction 

Recently, ultracold fermions have attracted much interest. One of the most interesting topics 
q | in the field is the superfluid state [Tj, where Cooper pairs composed of distinct components 

condense at low temperatures. Owing to the high controllability of the interaction strength 
and the particle number, interesting phenomena have been observed such as the BCS-BEC 
£> ' crossover [21 S] and the superfluid state in imbalanced systems [6]. These observations 

stimulate further experimental and theoretical investigations on fermionic systems. 

For optical lattice systems, which are formed by loading ultracold fermions into a periodic 
potential, various types of superfluid states have been proposed El [9]. One of the possible 

r~^! . candidates is the polarized superfluid state, which is naively expected to be realized in higher 

dimensional systems with imbalanced populations [10]. In our previous paper [11], we have 
investigated the attractive Hubbard model in infinite spatial dimensions to clarify how the 
polarized superfluid state is realized at low temperatures. However, dynamical properties were 
not yet studied in detail in the intermediate correlation and temperature region. In particular, 

S^ the question whether or not a pseudogap structure appears in the density of states above the 

critical temperature, has not been answered. To clarify this point, we study the formation of 
the superfluid gap in the system with and without spin- imbalanced populations. 

The paper is organized as follows. In §2, we introduce the model Hamiltonian and briefly 
summarize the DMFT framework. In §3, we study the attractive Hubbard model at quarter 
filling. By combining dynamical mean-field theory (DMFT) [X2^ [13| [14"1 [15] with the continuous 
time quantum Monte Carlo (CTQMC) method [16], we discuss how the polarized superfluid 
state is realized and the superfluid gap appears in the density of states at low temperatures. A 
brief summary is given in §4. 



2. Model and Method 

We consider a correlated fermion system with attractive interactions, which may be described 
by the Hubbard Hamiltonian, 

H= Yl [-*-(/" + ha) 5ij] c\ a c ja -U^ n^n^, (1) 
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where c%a (c ia ) is an annihilation (creation) operator of a fermion on the ith site with spin a, 
and riia = c ia Ci a . U is the onsite attractive interaction, t is the transfer integral between sites, 
/i is the chemical potential, and h is the magnetic field, which controls the spin imbalanced 
populations in the system. 

To study the attractive Hubbard model in infinite spatial dimensions [101 [TT1 [ITl HBl [191 [2D] 5 we 
make use of DMFT p2J[13l[TIl[T5]. In DMFT, the original lattice model is mapped to an effective 
impurity model, which accurately takes into account local particle correlations. The lattice 
Green's function is obtained via a self-consistency condition imposed on the impurity problem. 
This treatment is formally exact in infinite spatial dimensions, and DMFT has successfully been 
applied to strongly correlated fermion systems. When the superfluid state is directly treated 
in the framework of DMFT, the Green's function should be described by a 2 x 2 matrix. The 
self-consistency condition for the model [21] is given as, 

(D\ 2 
G l (icj n ) = (ioo n + h) a + fxa z - f — J a z G(iu n )a z , (2) 

where Go(iu) n )[G(iu n )] is the non-interacting (interacting) Green's function for the impurity 
model, w„,[= (2n + l)vrT] is the Matsubara frequency, T is the temperature, <5"o is the identity 
matrix, and a z is the z-component of the Pauli matrix. Here, we have used the semi-circular 
density of states, po(x) = 2/ttD^I — (x/D) 2 , where D is the half bandwidth. In the following, 
we set D as a unit of the energy. 

When low energy properties in strongly correlated systems are discussed in the framework of 
DMFT, an impurity solver is necessary to obtain the Green's function for the effective impurity 
model. Here we use the recently developed CTQMC technique, which has successfully been 
applied to the Hubbard model [11| I2"2"j 123] , the periodic Anderson model [21] , the Kondo lattice 
model |25j and the Holstein-Hubbard model [26J. In this paper, by using a CTQMC method in 
the continuous-time auxiliary field formulation [27] and extended to the Nambu formalism as 
an impurity solver, we discuss how the polarized superfluid state is realized in the system. The 
details of the method have been explained in Ref. |llj . 

3. Results 

We study the superfluid state in the Hubbard model at quarter filling to discuss how the gap 
structure appears in the density of states. A quarter filled system is chosen to study the attractive 
Hubbard model at a generic filling. The pair potential A(= (ci-j-cy.)) and the density of states 
Pa{w) in the case U = 2 are shown in Fig. [TJ When h = 0.0, as temperature is decreased, the pair 
potential is induced at T c ~ 0.095, as shown in Fig. [1] (a). This implies that a phase transition 
occurs from the normal metallic state to the superfluid state. A superfluid gap characteristic 
of this state clearly appears in the density of states at low temperatures, as shown in Fig. Q] 
(b). When the magnetic field is switched on, pairing correlations are suppressed and the normal 
metallic state becomes stable. In fact, the superfluid order parameter is decreased and the 
critical temperature is shifted to T c ~ 0.077 when h = 0.3, as shown in Fig. [T] (a). When a 
large magnetic field is applied, the polarized superfluid state is no longer stable even at very 
low temperatures. This instability has been discussed in Ref. [TT]. The density of states for 
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Figure 1. (a) Pair potential as a function of the temperature when h = 0.0 and h = 0.3. Solid 
and dashed lines represent the spectral functions for fermions with up (down) spin when h = 0.0 
(b) and h = 0.3 (c). 

each spin in the h = 0.3 case is shown in Fig. Q] (c). It is found that the density of states 
around the Fermi level decreases, with decrease of temperature from T c , similar to the case with 
h = 0.0. These results suggest that the superfluid gap appears only at T < T c . Therefore, we 
conclude that no pseudogap behavior is found at T > T c , at least, in the quarter filled Hubbard 
model with the inter ation U = 2. However, there still exists an open question whether or not 
a pseudogap behavior appears in other parameter regions such as the weak or strong coupling 
region and other fillings. Further systematic calculations are necessary to clarify this point, and 
will be discussed elsewhere. 

4. Summary 

We have investigated the attractive Hubbard model in infinite spatial dimensions, by combining 
dynamical mean-field theory with continuous-time quantum Monte Carlo simulations based on 
the Nambu formalism. We have calculated the superfluid order parameter and the density of 
states to discuss how the polarized superfluid state is realized at low temperatures. It was found 
that no pseudogap behavior appears above the critical temperature. 
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